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ABSTRACT 

We present a method for measuring the physical parameters of the coldest T-type brown dwarfs 
using low resolution near infrared spectra. By comparing H2O- and ^-sensitive spectral ratios be- 
tween empirical data and theoretical atmosphere models, and calibrating these ratios to measurements 
for the well-characterized 2-5 Gyr companion brown dwarf Gliese 570D, we derive estimates of the 
effective temperatures and surface gravities for 13 mid- and late-type field T dwarfs. We also deduce 
the first quantitative estimate of subsolar metallicity for the peculiar T dwarf 2MASS 0937+2931. 
Derived temperatures are consistent with prior estimates based on parallax and bolometric luminos- 
ity measurements, and examination of possible systematic effects indicate that the results are robust. 
Two recently discovered late-type T dwarfs, 2MASS 0939-2448 and 2MASS 1114-2618, both appear 
to be >50 K cooler than the latest-type T dwarf, 2MASS 0415-0935, and are potentially the coldest 
and least luminous brown dwarfs currently known. We find that, in general, higher surface gravity T 
dwarfs have lower effective temperatures and luminosities for a given spectral type, explaining pre- 
viously observed scatter in the T e ///spectral type relation for these objects. Masses, radii and ages 
are estimated for the T dwarfs in our sample using the evolutionary models of Burrows et al.; we also 
determine masses and radii independently for eight T dwarfs with measured luminosities. These two 
determinations are largely consistent, lending support to the validity of evolutionary models at late 
ages. Our method is well suited to large samples of faint brown dwarfs, and can ultimately be used 
to directly measure the substellar mass function and formation history in the Galaxy. 

Subject headings: stars: low mass, brown dwarfs - stars: fundamental parameters - 
stars: individual (2MASS J09373487+2931409, 2MASS J09393548-2448279, 
2MASS J11145133-2618235, Gliese 570D) 



1. INTRODUCTION 

The spectral energy distributions of the coldest 
known stars and brown dwarfs, L dwarfs and T 
dwarfs iKirkpatrick et all 119991: IBurgasser et~all l2002bt 
IGeballe et alJl2002j) . are complex, dominated by broad, 
overlapping gaseous and condensate molecular absorp- 
tion features. The strengths of these features de- 
pend on a combination of photospheric temperature, 
gas pressure and composition (e.g., Burrows & Sharp 
1999; Lodders & Fegley 2002), which in turn are re- 
lated to the effective temperature (T e //), surface grav- 
ity (g) and metallicity ([M/H]) of a brown dwarf. 
Nonequilibrium effects (e.g., vertical mixing, cloud cov- 
erag e) may also play an i mportant role in molecu- 
lar llFeglev fc Loddersl 119961 iSaumon et alJ 120031) and 
conden sate (jAckermarifeMarle^l27)0^^ et alJ 

l2002aft abundances. The combined influence of these pa- 

1 Visiting Astronomer at the Infrared Telescope Facility, which is 
operated by the University of Hawaii under Cooperative Agreement 
NCC 5-538 with the National Aeronautics and Space Administra- 
tion, Office of Space Science, Planetary Astronomy Program. 



ramcters on the spectra of L and T dwarfs is only begin- 
ning to be explored through the study of low mass, sub- 
stellar objects in young clusters and st ellar associations 
llLucas et aTll200lt iGorlova et alJl2003l: iMcGovern et all 
20041) and ultracool subdwarfs fBurgasseret^dTT2003at 
Scho lz. Lodieu fc McCaughrearl2004[) . although svstem- 
atic studies have yet to be achieved. 

Disentangling the properties of T e ff, surface grav- 
ity and metallicity is a principal goal of substellar as- 
trophysics. These parameters can be used to infer 
masses, radii and ages for individual sources (e.g., Mo- 
hanty, Jaywardhana & Basri 2004), allowing, in the long 
term, direct measurement of the substellar mass func- 
tion (MF) and star formatio n history for field objects 
in the Solar Neighborh ood ijChabrierl 120031 iBurgasseil 
120041 lAllen et all 12005ft . In the short term, T e// , g 
and [M/H] measurements for young cluster or compan- 
ion brown dwarfs ena ble tests of evolutionary models 
l)Mohantv et alJl2004aft . 

Gravity and metallicity effects are particularly rele- 
vant for interpreting the spectral energy distributions of 
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the coldest T dwarfs, spectral t ypes T6 and later. These 
objects, with T e// < 1000 K jGolimowski et al J 120041 ). 
lack the complicating influence of photospheric conden- 
sates common in late-typ e M dwarfs, L dwarfs and the 
earliest-type T dwarfs (iTsuii. Ohnaka. fc Aokil 119961 
lAc.kerman fc Marlevl 120011 lAlla.rH et a, U I2001D . They 
exhibit good c orrela ti on between spectral type and T P f f 
llDahn et al.l l2002t iTinnev. Bu rgasser. fc Kirk patricld 
20031 IGolimowski et alJ 120041 INakaiima et all 120041 



our method for various brown dwarf studies. Results are 
summarized in § 7. 



VrbaetldTl27)o"4ji . Surface gravity and metallicity effects 



are therefore readily distinguished by the presence 
of spectral or photometric anomalies. One case in 
point is the pe culiar T6 dwarf 2MASS 0937+2931 2 
ijBurgasser et alJ l2002b|) . a brown dwarf believed to 
have a high surface gray i ty and/ or subsolar met allicity 
(IBurgasser et alJ I2002bl I2003bl IBurrows et all 120021: 
IKnapp et al] 120041) . 2MASS 0937+2931 is 0.5-1.0 mag 
bluer than similarly-classified T dwarfs; and its spectrum 
exhibits a suppressed if-band peak, an extremely red 
0.8-1.0 /im spectral slope, enhanced FeH absorption at 
0.99 fim, and an absence of K I doublet lines at 1.17 and 
1.25 /im, all unusual for a mid- type T dwarf. Several 
other late-type T dwarfs exhibit similar color and spec- 
tral peculiarities (iBu rgassen McElwain. fc K irkpatrickl 
120031: iBTugasser et alJl2004bt IKnapp et al.ll2004fl . How- 



ever, quantitative analysis of these effects, in the form 
of specific surface gr avity and meta l licity measurements, 
has been limited ijBurrows et al.l l2002t IKnapp et al.l 
l200l . 

In this article, we present a method for disentangling 
T e ff, 9 and [M/H] effects in the near infrared spec- 
tra of the latest-type T dwarfs. Our method, based 
on the comparison of calibrated near infrared flux ra- 
tios measured on low resolution spectral data and theo- 
retical models, yields strong constraints on these physi- 
cal parameters and a means of estimating masses, radii 
and ages for individual field brown dwarfs. In § 2 we 
describe the sample and spectrosco pic observations ob- 
tained with the SpeX spectrograph (jRavner et alJl200"3|) 
mounted on the 3m NASA Infrared Telescope Facility 
(IRTF) . We identify and compare spectral variations ob- 
served in these low resolution near infrared spectra, and 
discuss qualitatively how these features are associated 
with differences in T e //, g and [M/H]. In § 3, we examine 
these same effects with theoretical models, and charac- 
terize spectral trends. In § 4, we describe our method, 
and present T e ff and log<? estimates for 13 field brown 
dwarfs and constraints for two others; we also deduce 
subsolar metallicity estimates for two sources including 
2MASS 0937+2931. We demonstrate the consistency of 
our T e ff values with previous determinations based on 
parallax and luminosity measurements, and examine po- 
tential systematic effects. In § 5, we derive mass, radius 
and age estim ates for our T dwarfs using the evolution- 
ary models of IBurrows et al.l l|1997j) : and independently 
determine masses and radii for eight sources with pub- 
lished luminosity measurements. We discuss the results 
in § 6, focusing on new insights on the T e /j/spectral 
type relation for T dwarfs and potential applications of 

2 Source designations in this article are abbreviated in the man- 
ner 2MASS hhmmiddmm; the suffix is the sexagesimal Right As- 
cension (hours and minutes) and declination (degrees and arcmin- 
utes) at J2000 equinox. Full designations are provided in Table HI 



2. OBSERVATIONS 
2.1. The Sample 
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the Two Micron All Sky Survey l|Cutri et all 12001 
hereafter 2 MASS) and the Sloan Digital Sky Survey 
(|York et all 120001 hereafter SDSS). The empirical 
properties of these sources are listed in Table ^ The 
sample was selected to span types T5 .5 to T8, based 
on the unified classification scheme of Burgasscretal 



(2005), and excludes known binaries IjBurgasser et al 
2003c, Burgasser et al. in prep.). Eig ht of these 
object s have parallax measurements fro m Dah n et alJ 
(2002 ): ITinnev. Burgasser. fc Kirkpatrick! l|2003l) : and 
Vrb a et all p004): all but one has a reported proper 
motion. Apparent 2MASS J-band magnitudes for these 
sources range from 14.7 to 16.3 mag. 

2.2. Near Infrared Spectroscopy 

Six of the T dwarfs in our sample 
2MASS 0034+0523, 2MASS 0243-2453, 2MASS 0415- 
0935, 2MASS 1231+0847, Gliese 570D and 2MASS 2228- 
4310 - have been pr eviously observed with SpeX 
(Burgasse r et al.l l2004b). The remaining sources were 
observed during three runs on 2004 March 11-12, 2004 
July 23 and 2004 September 7 (UT). A log of obser- 
vations is provided in Table El Conditions during the 
March run were clear and dry with typical seeing of 0'.'7. 
Conditions during July were also clear with excellent 
seeing (0'.'4-0'.'7). Light cirrus was present during the 
September observations, but seeing was again excellent 
(0 / /5-Q'/7). 

Spectral data for all of the sources in our sample (in- 
cluding those previously observed) were obtained using 
the SpeX prism dispersed mode, which provides low res- 
olution 0.7-2.5 fxm spectra in a single order. This set- 
ting minimizes spectral color errors commonly incurred 
through order stitching (e.g., McLean et al. 2003), yield- 
ing an accurate measure of the broad band spectral en- 
ergy distribution. For all observations, the 0'.'5 slit was 
employed and rotated to the parallactic angle, result- 
ing in a spectral resolution A/AA « 150 and dispersion 
across the chip of 20-30 A pixel^ 1 . Multiple exposures 
of 180 s were obtained in an ABBA dither pattern along 
the slit. Flux calibration was made through observa- 
tions of nearby AO V stars obtained immediately be- 
fore or after the target observation and at similar air- 
masses (Asecz < 0.1). Internal flat field and Ar arc 
lamps were observed after each flux calibrator star for 
pixel response and wavelength calibration. All spec- 
tral data were reduced using the SpeXtool package ver- 
sion 3.2 l|Vacca et al J l2003t iCushing. Vacca. fc Ravnerl 
2004) using standard settings. Further details on the 
experimental design an d data reduction are given in 
IBurgasser et alJ l)2004hl) . 

The reduced spectra of the newly observed T dwarfs 
are shown in Figure ^ Readily apparent are the deep 
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molecular bands of H2O and CH4 that shape the 1.05 
(Y-band), 1.27 (J-band), 1.6 (ii-band) and 2.1 /im 
(if -band) flux peaks, the defining features of T dwarf 
near infrared spectra. The spectra are also shaped 
by the pressure-broadened red wings of the 0.77 //m 
K I doublet shortward of 1 /im and collision-induced 
H2 absorption at if -band, both of which are discussed 
in detail below. Finer atomic line features, including 
the 1.17 and 1.25 /im K I doublets, are unresolved in 
these data. Further discussion on the spectral charac- 



teristics of T d w arfs can be fo u nd in |Burgasser^^d 



2003H): ICeT^nil EQOai: IMcLear, et 
Kna pp et all ( 2004); Nakaii maeTaLl ( 200 



al 



Cushing, Ra vner. fc Vaccal l)2005j) : and [Kirkpatric 



2.3. Spectral Signatures of Surface Gravity and 
Metallicity 

Variations in the near infrared spectral features of T 
dwarfs are generally synchronized with spectral type - 
later subtypes exhibit both stronger H2O and CH4 bands 
and bluer near infrared colors. However, slight devi- 
ations to these trends exist, and are apparent when 
one compares sources with similar spectral types, as 
in Figure |21 Displayed in the left panel of this fig- 
ure are the normalized spectra of three T6/T6.5 dwarfs, 
2MASS 0937+2931, SDSS 1346-0031 and 2MASS 
2228—4310, overlai n on that of the T6 spe ctral standard 
SDSS 1624+0029 ijBureasser et all 12005]) . While H 2 
and CH4 bands are similar among these spectra, clear 
differences are seen in the relative brightness of the K- 
band flux peak and the shape of the Y"-band peak. In 
particular, 2MASS 0937+2931 exhibits weaker if-band 
emission and a broader Y"-band flux peak as compared to 
SDSS 1624+0029, while 2MASS 2228-4310 has stronger 
if-band emission. Similar deviations are also seen 
among the three T7.5/T8 dwarfs 2MASS 0939-2448, 
2MASS 1114-2618 and 2MASS 1217-0311 when com- 
pared to the similarly classified Gliese 570D. 

What gives rise to these deviations? Shortward of 
the y-band spectral peak, the dominant absorbers in T 
dwarf spectra are the pressure-broadened wings of the 
K I and Na I fundamental doublet lines centered at 
0.77 a nd 0.59 
1999f jBurrgw^ 
120031 Burrows fc Volobuvevl |2003j). These features 
strengthen with late r spectral type throughout the L and 
T dw arf sequences (iKirkpatrick et alJ 119991: iReid et aP 
120001 iBurgasser et al.ll2003b|) . The broad wings of the 
alkali lines, induced by kinematic perturbations by other 
chemical species (most importantly H2 and He), are en- 
hanced in higher pressure (P) and higher density atmo- 
spheres. As atmospheric pressure scales as dP/dr ~ 
P/t oc g/nR (where r is the optical depth and kr the 
Rosseland mean opacity), higher pressure photospheres 
(r = 2/3) are achieved in brown dwarfs with higher 
surface gravities and/or metal-deficient atmospheres (re- 
duced kr). For these sources, line broadenin g theory 
llAllard et alJl2003l lB urrows fc Volobuvevl2003j) predicts 
the strongest absorption close to the line centers, result- 
ing in steep 0.8-1.0 /im spectral slopes due to red wing 
of the K I doublet. 

The if -band peak, while molded by H2O and 
CH4 bands at 1.8 and 2.2 /im, is dominated 



/im, respectively llTsuii. Qhnaka. fc Aoki 
vs. Marlev. fc Shard 1200(1 lAllard et al. 



by another press ure-sens it ive fe at ure, collision-i n duced 
H 2 absorption lILinskvl 119691: ISaumon et"aH 119941 
Borvsow, Jor gensen. fc Zhend 119971) . The ind~ced 1-0 
quadrupolar moment of this molecule produces a broad, 
featureless absorption centered near 2.1 /im. Like the 
K I wings, H2 absorption arises from kinematic pertur- 
bations and is therefore enhanced in the higher pressure 
and higher density atmospheres present on high surface 
gravity and/or low metallicity brown dwarfs. 

While deviations in the strengths of the K I and 
H2 features have previously been li nked to gravity 
and metallicity variations in T dwarfs iBurga sser et all 
2002bl I2003bt IBurrows et all 120021: ILeggett et alJ 120031: 



Knapp et aiT l2004). Figure [21 demonstrates that these 
features are correlated. The steeper K I wings 
and enhanced H2 absorption exhibited in the spec- 
tra of 2MASS 0937+2931, 2MASS 0939-2448 and 
2MASS 1114-2618 are both indicative of higher pres- 
sure photospheres; while the weaker H2 absorption in the 
spectra of 2MASS 2228-4310 and 2MASS 1217-0311 in- 
dicate low pressure photospheres. 

In contrast, the congruence of the CH4 and H2O bands 
for similarly classified T dwarfs suggests that gravity and 
metallicity effects for these features are minimal. The ob- 
served correlation between T e // and spectral type, the 
latter based on the strengths of the molecular bands, 
links H2O and CH4 to temperature. However, gas pres- 
sure does regulate the atmospheric abundance of CH4 
and H 2O in the principle reac t ion CO + 3H2 CH4 + 
H 2 (lFeglev fc LoddersllT996l IBurrows fc Shard 1 1999ft . 
while metallicity modulates bo th CH4 and H2O abun- 
dances (jLodders fc Feglevl2002|) . Hence, nearly all of the 
major absorption features in T dwarf spectra are affected 
in some manner by T e ff, 9 and [M/H]. 

3. SPECTRAL MODELS 

To further investigate the physical origins of the 
spectral peculiarities described above, we have ex- 
amined a new suite of brown dwarf spectral mod- 
els that incorporate differences in T e ff, surface grav- 
ity and metallicity. The models, develop e d by 
the Tucson group llBurrows. Marlev. fc Shard 200C ; 
Burrows et al.ll2002t BuT ows. Sudars kv fc Lunind l2003: 



Burrows. Sudarskv fc Hubenyjl2005|) . are self-consistent , 
non-gray atmospheres incorpora ting up-to-date molec- 
ular opacities as described in IBurrows et al.l (2001). 
The atmospheres are assumed to be free of conden- 
sate dust species, consistent with prior modellin g results 
l)Tsuii. Qhnaka. fc Aokill999HAllard et alJ2001|) . follow- 
ing the prescription of co ndensate rainout as described in 
IBurrows fc Shard (|1999l see also Lodders 1999). Modi- 
fied Lorentzian profiles with an ad-hoc, smooth cutoff arc 
used to model the line broadening of the Na I and K I 



doublets l)Burrows. Marlev. fc Sharoll200l . Nonequilib 



— - - — ft) |^ (mmj, - — . 

rium mixing effects llSaumon et al.l l2003j) are not con- 



sidered. A ful l des cription of these models is given in 
IBurrows et aD l|2002ft . 

Figure compares three sequences of these spectral 
models varying T e //, g and [M/H], respectively. The 
spectral resolution of the models has been degraded using 
a Gaussian kernel to match the resolution of the SpeX 
prism data. There is general agreement in the overall 
spectral morphologies of the models and observed data; 
however, important discrepancies are present. Most 
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prominent of these is the shape of the 1.6 /jm CH4 band, 
reflecting continued deficiencies in the near infrared opac- 
ities of this molecule for which only low temperature 
(300 K) laboratory measurements have been obtained 
(Saumon et al. 2000; however, see Homeier et al. 2003). 
These opacities also detrimentally affect absorption fea- 
tures at 1.1 and 1.3 /urn, although bands at 2.2 and 3.3 /im 
have been found to be adequately reproduced (M. Cush- 
ing, 2005, priv. comm.). In addition, the line-broadening 
theory employ ed by these models predates more recent 
calculations bv lBurrows fc Volobuvevl l)2003|) . which pre- 
dict a sharper cutoff for the red K I wing at 1.0 /jm, 
in contrast to the modified Lorentzian profile used here 
which is relatively flat over the 0.9-1.0 /im waveband. 
As a result, alkali opacity at shorter wavelengths in the 
models is reduced (note the stronger 0.92 /im H2O as 
compared to the data) while the 1.05 (im F-band peak 
is more suppressed in the models than observed. 

The top panel of Figure[3]shows temperature variations 
for solar metallicity and log g = 5.0 cm s~ 2 models. The 
trends in this sequence reflect those observed in T dwarf 
spectra as a function of spectral type; i.e., strengthening 
H2O and CH4 bands producing more acute triangular 
flux peaks, and stronger absorption shortward of 1 /im 
and at A"-band. The increasing depths of the molecular 
bands with cooler effective temperatures is largely a con- 
sequence of the increased column depth, and therefore 
total opacity, of the associated gas species. 

The middle panel compares surface gravity variations 
in the T e ff = 800 K, solar metallicity models. Here, 
spectral variations are strongest at the Y- and A-band 
flux peaks, although H2O and CH4 absorptions at 1.1 and 
1.3 /im, and the iAband peak, are also affected. The A- 
band peak is suppressed in higher surface gravity models, 
consistent with the enhanced H2 absorption expected in 
higher pressure photospheres. On the other hand, sur- 
face gravity variations at the V-band peak do not match 
the observed trends. The highest surface gravity model 
exhibits reduced 1.05 /im flux relative to 1.27 /im, con- 
trary to the brighter and broadened K-band peaks seen 
in the empirical data. We attribute this discrepancy to 
the outdated alkali line broadening profile employed in 
these models, and leave analysis of this feature to a fu- 
ture study. 

The similarity of surface gravity modulations of the 
A-band flux peak in the theoretical models to variations 
in the spectral data shown in Figure □ is highlighted 111 
Figure 01 which shows a similar sequence of T6/T6.5 
and T7.5/T8 dwarfs but overlain on the T e// = 1000 
and 800 K solar metallicity models, respectively. In- 
equities in the 1.6 /xm CH4 band notwithstanding, the 
spectral models reproduce reasonably well the relative 
variations observed in the A-band flux peaks, although 
both 2MASS 0937+2931 and 2MASS 0939-2448 exhibit 
stronger A-band suppression than the highest surface 
gravities permitted by the models. 

The bottom panel of Figure |3 compares metallicity 
effects at fixed surface gravity (log g = 5.5) and tem- 
perature (T eff = 800 K) for [M/H] = 0, -0.5 and -1. 
The higher surface gravity examined here is appropriate 
for old brown dwarf members of the metal-poor Galac- 
tic thick disk and halo populations. Spectral variations 
are far more extreme in this case. The A-band peak 
is suppressed at lower mctallicitics, as expected for en- 



hanced H2 absorption. At the same time, emergent flux 
appears to be relatively enhanced shortward of 1 /im, due 
to reduced Na and K abundances and their correspond- 
ing opacities. The specific shape of the spectrum at these 
wavelengths should be treated with caution, however, 
given the older line broadening theory used in the models. 
Shifts in the J- and iJ-band peak wavelengths are due to 
reduced H2O abundances and increased collision- induced 
H2 absorption extending into the ff-band. Overall, these 
spectral variations are more substantial than those seen 
in the empirical data, a sign that significant subsolar 
metallicities ([M/H] < —0.5) are not present among 
the T dwarfs examined here. However, the broadened 
K-band peaks and strong A-band suppression in the 
spectra of 2MASS 0937+2931, 2MASS 0939-2448 and 
2MASS 1114-2618 do hint at slightly subsolar metallic- 
ities for these sources. 

4. MEASURING PHYSICAL PROPERTIES FROM T DWARF 
SPECTRA 

4.1. The Method 

The spectral models confirm that the pressure-sensitive 
H2 and K I absorptions are more strongly influenced by 
surface gravity effects than the absorption bands of H 2 
and CH4, while the latter vary more strongly with T e ff, 
at least for the temperature regime considered here. By 
contrasting the strengths of these features, it should be 
possible in principle to extract the effective temperatures 
and surface gravities of these objects. In practice, this 
pursuit has proven problematic due to persistent inade- 
quacies in molecular opacities and corresponding system- 
atic errors in the models (e.g., Burgasser et al. 2004a). 
What is required is a means of calibrating the spectral 
models using one or more empirical fiducials. 

Fort unately, such a fiducial exists in the T dwarf Gliesc 
570D (Burgass er et al.ll2000jh This widely-separated (p 
— 258" = 1530 AU), common proper motion brown dwarf 
companion to the nearby (5.91+0.06 pc; ESA 1997) 
Gliese 570ABC system has both distance and luminos- 
ity measurements that are empirically well constrained. 
The age of this system is estimated to be 2-5 Gyr based 
on a comparison of age, activity and ki nematic relations 
for t h e K and M stellar components l)Burgasser et alJ 
l2000t lOeballe et aDl200l . The K4 V primary has a 
near-solar metallicity ([Fe/H] = 0.00-0.16; Feltzing k 
Gustafsson 1998; Thoren & Feltzing 2000; Allende Pri- 
eto & Lambert 2000). Assuming coevality, Gliese 570D 
is one of the fe w T dwarfs with both age and metallic- 
ity constraints. iGeballe et alJ l)2001|) derive fairly pre- 
cise temperature (T e / / = 784-824 K) and surface grav- 
ity (g = (1-2) xlO 5 cm s~ 2 ) estimates, the former based 
on the measured luminosity and a model-dependent ra- 
diu s, the latter based on brown dwarf evolutionary mod- 
els ijBurrowset al .111 997ft. 

Using Gliese 570D as our empirical fiducial, our pro- 
cedure was then as follows. The strengths of the major 
H2O bands and relative fluxes of the spectral peaks were 
measured for both the empirical data sample and theo- 
retical models using the following ratios: 



/ Fl, 26-1.285 ' 

H 2 0-H= {^ 48 - 1 - 62 , (2) 

J -fl. 56-1. 60 
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Y/J 



005-1.045 



J*. 



25-1.29 



and 



v , j J F-2. 06-2. 10 

K/J= T - 

J #1. 25-1.29 

K/H= { F J 06 - 210 , 

J -Fl. 56-1. 60 



(3) 



(4) 



(5) 



where J Fx 1 -\ 2 denotes the integrated flux between 
wavelengths Ai and A2 (in microns). The first two ra- 
tios are identical to those defined for the near infrared 
classification of T dwarfs; the s pectral region sampled 
by these is shown in Figure 5 of Burgas ser et al.l (2005, 
note that ratios sampling the poorly modelled CH4 bands 
are not considered here). The K/J ratio has also been 
used previously to examine v ariations in H2 absorption 
l|Burgasser et al.ll2004bll2005|) . The Y/J and K/H color 
ratios are defined here for first time. Measurements of 
these ratios on the empirical data are given in Table |3| 

A series of H2O— J and K/H ratios for the solar metal- 
licity models are shown in Figure Both ratios vary ac- 
cording to differences in T e ff and log g, although gravity 
variations are stronger in the K/H color ratio. Gravity 
variations also affect the two ratios in opposite ways; 
higher gravity models yield large H 2 0— J values (imply- 
ing weaker absorption) and smaller K/H values (imply- 
ing weaker if-band emission). Similar trends are seen 
with the H2O— H and K/J ratios, respectively. All of 
these trends are consistent with the qualitative proper- 
ties of the model spectra shown in Figure |31 

Calibration of the model ratios was achieved by cor- 
recting these values to those measured from the SpeX 
spectrum of Gliese 570D. Adopting T e ff — 800 K and 
logg = 5.1 for this source, we computed the correspond- 
ing model ratios by linear interpolation. Correction fac- 
tors, listed at the bottom of Table |3 were defined as 
the ratio of the spectral data measurement to the model 
value. For four of the ratios, model values differ by less 
than 20%; the H2O— H index, on the other hand, requires 
a 60% correction. The correction factors were applied to 
all of the solar metallicity model ratios, regardless of T e f / 
or gravity, and therefore represents a first order calibra- 
tion of the models. 

In principle, any of the ratios defined in Eqns. 1-5 
could be used for comparison to the spectral data. We re- 
strict our primary analysis to the H 2 0— J index, which 
requires a smaller correction factor than H2O— H; and 
the K/H index, which gives a quantitative measure of 
the behavior demonstrated in Figure 01 The corrected 
model ratios were resampled in steps of 20 K in T e ff 
and 0.1 dex in logg by linear interpolation. Then, for 
each spectrum, we identified the region in T e f / and log g 
phase space for which the corrected model ratios agreed 
with the empirical ratios, assuming a 10% uncertainty 
(see § 4.3). 

Figures El and illustrate these matches. The H2O— J 
and K/H ratios each constrain a set of T e ff and logg 
values that span the model parameter space diagonally; 
e.g., agreement in the H 2 0— J ratios span low tempera- 
tures and high surface gravities to high temperatures and 
low surface gravities. The K/H ratios match an orthog- 
onal phase space. The intersection of these phase spaces 



provides an unambiguous constraint on both T e ff and 
logg- 

4.2. Results 
4.2.1. Teff and g Estimates 

Table 0] lists the ranges of T e f f and log g constrained 
by the two ratios for each source in our sample. For 
13 of the 16 sources (including Gliese 570D), these val- 
ues are well defined, with typical uncertainties of 40- 
60 K in T e ff and 0.1-0.3 dex in logg. For three of 
the T dwarfs, 2MASS 0937+2931, 2MASS 0939-2448 
and 2MASS 1114—2618, no phase space intersection was 
found. The case of 2MASS 0937+2931 is discussed 
in further detail below. For 2MASS 0939-2448 and 
2MASS 1114-2618, close examination of FiguresEJandQ 
suggests that phase space intersections are possible at 
lower T e //s than those spanned by our model set; i.e., 
for T e ff < 700 K. This is intriguing, since we derive 
T e// = 740-760 K for the T8 2MASS 0415-0935, cur- 
rently t he coldest and lowest luminosity brown dwarf 
known iBurgasser et all 120024 IGolimowski et alJ 120041 
iVrba et aLll2004f) . Parallax measurements can determine 
whether 2MASS 0939-2448 and 2MASS 1114-2618 are 
in fact colder and fainter brown dwarfs. 

Are these temperatures and surface gravities reason- 
able? Eight of the T dw arfs in our sample have prior 
T F jf determinati ons from IGolimowski et alJ (|2004f) and 
IVrba et all l)2004l) based on parallax and bolometric lu- 
minosity (L) measurements; these values are listed in 
Tabled In all eight cases, our derived T e //s are consis- 
tent. This agreement may have much to do wi t h the large 
T e _ f f estim a te ran ges from IGolimowski et alJ i|2004|) and 
IVrba et aP l)2004|) . as high as 300 K, due to uncertain- 
ties in the radii adopted to compute T e ff from L. Our 
T e f f estimates are typically in the middle or high end of 
the ranges from these studies. The only surface gravity 
es timates for fie l d T d warfs reported to date are those 
of iKnanp et al.1 l|2004fl . based on a comp arison of near 
infrare d colors to atmosphere models by iMarlev et alJ 
(2002). In this case, we find that our estimates are sys- 
tematically 0.3-0.5 dex higher than the Knapp et al. val- 
ues. As the latter are stated without uncertainties, we 
simply point out this discrepancy for further study. 

4.2.2. A Subsolar Metallicity for 2MASS 0937+2931 

The parameter spaces for 2MASS 0937+2931 do not 
intersect in Figure H but it appears that they would if 
higher surface gravities were modelled. However, sur- 
face gravities larger than logq = 5.5 ar e restricted by 
the interior physics l|Burrows et al.lll9 971 An alternate 
hypothesis is that the spectrum is influenced by a third 
parameter, namely metallicity. 

We can quantitatively test this case by introducing 
metallicity variations into the model set. Applying the 
same corrections to the model ratios as above (i.e., as- 
suming Gliese 570D has [M/H] = 0), linearly interpolat- 
ing between the [M/H] = and [M/H] = -0.5 models 
in 0.1 dex steps for 700 < T eff < 1200 K and 5.0 < 
logg < 5.5, and performing the same comparative anal- 
ysis, we derive the series of parameter phase spaces (ef- 
fectively slices of a three dimensional parameter phase 
volume) shown in Figure [SJ We find that, for the case of 
2MASS 0937+2931, the phase spaces intersect for metal- 
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licities -0.1 < [M/H] < -0.4, 780 < T e// < 860 and 5.0 

< log g < 5.5. Lower metallicities may also be feasible at 
lower surface gravities. One other source in our sample, 
2MASS 0034+0523, also exhibits intersecting parameter 
spaces for slightly subsolar metallicities, —0.2 < [M/H] 

< 0. This object has also been noted for its strong K- 
band suppression l)Burgasser et al.ll2004b|) . 

Our analysis represents the first quantitative con- 
straints of metallicity for a brown dwarf, and are con- 
sistent with prior qualitative conclusions. However, the 
derived values should be treated with caution for rea- 
sons other than the fidelity of the theoretical models. 
Constraining the three parameters T e ff, 9 and [M/H] 
cannot be done unambiguously using only two spectral 
ratios; at least one additional constraint is required. A 
promising candidate is the Y / J ratio, as the metal-poor 
models shown in Figure^exhibit significant variations at 
these short wavelengths. We defer more thorough exam- 
ination of this third index until such time as the current 
generation of spectral models incorporate a more rigor- 
ous line broadening theory. 

4.3. Assessing Systematic Effects 

While the T e / / and log g estimates made here are rea- 
sonable and generally consistent with prior work, it is 
important to identify and characterize any source of sys- 
tematic error that may skew the results. Such system- 
atic effects can arise from the calibration or quality of 
the spectral data, the choice and calibration of the spec- 
tral ratios used and limitations of the spectral models 
themselves. We examine these effects here in detail. 

4.3.1. Flux Calibration of the Spectral Data 

Accurate measurement of color ratios such as K/H as- 
sumes that the spectral data portrays the true color of 
the source, which brings into question possible redden- 
ing of the observed spectrum and the accuracy of the 
overall flux calibration. Interstellar reddening can gen- 
erally be ruled out for our sample as all of the sources 
lie at distances of ~20 pc or less. A more local source 
of reddening, differential color refraction through our at- 
mosphere, has been mitigated by observing the sources 
with the slit aligned at the parallactic angle. We there- 
fore assume that both of these effects are negligible. 

Systematic errors incurred in the flux calibration can 
be quantified by comparing spectrophotometric colors 
for the data to published photometry. We examined 
J — H, H — K and J — K colors on the Mauna 
Kea Observatory system (MKO; Simons & Tokunaga 
2002; Tok unaga, Simons &: Vac ca 2002) using photom - 
etry from IGeballe et alJ l|2001fl : ILeeeett et all l|2002j) ; 
and IKnapp et all l|2004|) : and J - H, H - K s and 
J - K s colors from 2MASS (Table GJ. Spectrophoto- 
metric colors were determined by integrating the ap- 
propriate filter profile (combined with telescope and in- 
strumental optical response curves for 2MASS photom- 
etry; see § Ill.l.b.i in Cutri et al. 2003) over the near 
infrared spec tra of each T dwarf and that of the AO 
V star Vega ({Bergeron. Wesemael. fc Beauchanml 119951 
see also Stephens & Lcggctt 2003). We found no 
systematic differences for any of the photometric and 
spectrophotometric colors on both systems, and typi- 
cal deviations were 5% or less for the more accurate 
MKO photometry. For those few sources with color 



offsets significantly greater than their photometric un- 
certainties (2MASS 0243-2453, 2MASS 0727+1710 and 
SDSS 1758+4633 have 3a deviations in MKO J-K), dif- 
ferences were at most 15%. Hence, we conclude that the 
10% uncertainties adopted for the color ratios adequately 
compensates for uncertainties in the flux calibration. 

4.3.2. Spectral Noise 

Molecular band ratios are generally insensitive to color 
errors in the overall spectrum, but deep absorption bands 
can be affected by spectral noise. This is manifested 
by variations in the measured flux at the bottom of the 
bands where signal-to- noise (S/N) is minimal. To explore 
the impact of this effect, we performed a Monte Carlo ex- 
periment, measuring the H2O— J ratio on the combined 
spectrum of Gliese 570D plus a Gaussian noise compo- 
nent scaled to S/N = 10-200 at the J-band peak. A total 
of 10 000 trials were run over a uniform range of S/N. 
We found no systematic deviations in the H2O— J ratios 
for S/N > 10, but scatter among the values increased in 
the noisier spectra, approaching 10% for S/N = 40. All 
of our spectra have S/N > 50 at the J-band peak with 
the exception of SDSS 1110+0116, which has S/N « 20. 
Hence, we find that our adopted 10% uncertainties for 
the measured H2O ratios incorporate noise effects suf- 
ficiently, although derived values for SDSS 1110+0116 
may be more uncertain. 

4.3.3. Choice of Spectral ratios 

The use of the H2O— J and K/H ratios in our analysis 
above was justified by the sensitivity of these ratios to 
T e ff and g variations in the models, the magnitude of 
the calibration correction required and the fidelity of the 
models in these spectral regions. However, we can also 
consider how the results change if a different set of ratios 
are employed. We repeated our analysis with four pair- 
ings among the H2O— J, H2O— H, K/H and K/J ratios. 
Again, we found no significant or systematic differences 
in the derived T e j f and log g values amongst our sample, 
although parameters for individual sources differed by as 
much as 160 K and 0.6 dex, respectively. Typical devi- 
ations were of order 65 K in T e ff and 0.15 dex in logg, 
which we adopt as estimates of systematic uncertainty. 

4.3.4. Calibration of the Spectral ratios 

The calibration of the spectral ratios hinges largely 
on the assumed physical properties for our calibrator 
source Gliese 570D. But how sensitive are the results to 
these adopted parameters? By varying the assumed T e f f 
and log g for Gliese 570D by ±25 K and ±0.1 dex, re- 
spectively, consisten t with the range of values found by 
IGeballe et al] lf200T|) . we found mean offsets of ±35 K 
and ±0.1 dex in the derived parameters for the field 
sources. These offsets were independent; changing the 
adopted T e ff of Gliese 570D had no impact on the de- 
rived surface gravities, and vice versa. Hence, we es- 
timate that additional systematic uncertainties of 35 K 
and 0.1 dex in T e ff and log g, respectively, accommodate 
uncertainties in the physical parameters of our compari- 
son source. 

4.3.5. Choice of Spectral Models 
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Many of the potential systematic effects involving the 
choice and calibration of the spectral ratios would be 
eliminated if the models accurately reproduced the ob- 
served spectra. As limitations in the opacities prevent 
this, we must also consider how dependent our results 
are on the choice of spectral models used. We there- 
fore repeated our analysis using the COND models of 
lAllard et alJ 1)200 ID . Calibration of the spectral ratios 
was performed in the same manner as with the Tucson 
models, yielding somewhat different correction values. 
While deviations typically of order 50 K and 0.1 dex were 
found when comparing derived T e ff and log g between 
the model sets, these deviations were not systematic. 
Hence, we conclude that the choice of spectral model 
does not systematically change our results. 

In summary, we find no systematic deviations in our 
method that would lead to skewed estimates of T e // 
and log g, although systematic uncertainties of order 50- 
100 K and 0.1-0.25 dex may be present. 

5. MASS, RADIUS AND AGE ESTIMATES 
5.1. Evolutionary Models 

According to brown dwarf evolutionary theory, the 
temperature and surface gravity of a brown dwarf at a 
given time is directly related to its mass and age as- 
suming a given composition 3 . This implies that if T e / / 
and g are known, estimates for the latter, more funda- 
mental physical properties could be derived on an indi- 
vidual basis. This is demonstrated graphically in Fig- 
ure El which compares the derived T e / / and log g val- 
ues for 14 of our sources (excluding 2MASS 0939-2448 
and 2MASS 111 4-2618) to solar met allicity evolution- 
ary models from lBurrows et"al"l l)1997t) . Table El lists the 
range of masses and ages, and the corresponding radii, 
derived from this comparison. Our sample appears to 
span a broad range of masses (0.02-0.07 M Q ) and ages 
(<1-10 Gyr), consistent with a random sample drawn 
from the local Galactic environment. At late ages, brown 
dwarf radii are fairly constant, so mass and surface grav- 
ity are almost directly related. Hence, our lowest (high- 
est) surface gravity objects are also the least (most) mas- 
sive and youngest (oldest). Note in particular the place- 
ment of 2MASS 0937+2931, which appears to be the 
most massive and oldest in the sample. On the other 
hand, our analysis suggests that the low surface gravity 
T dwarfs SDSS 1758+4633 and 2MASS 2228-4310 may 
have masses less than 0.03 M Q and ages less than 1 Gyr. 

How reliable are these estimates? An independent 
check of the derived ages can be made by examining 
the kinematics of our sample. Three dimensional space 
velocities have not yet been measured for the T dwarfs 
examined here; a rough analysis can be made, however, 
by examining their tangential velocities, Vtan- Of the 14 
sources in our sample with age estimates, 13 have proper 
motion measurements and eight have parallax measure- 
ments. For 2MASS 0034+0523, 2MASS 0050-3322, 
2MASS 1231+0847 and 2MASS 2228-4310 we 
adopted spectrophotometric distance estimates 

3 We ignore for this discussion metallicity effects in the evolution 
of a brown dwarf, in addition to other variations related to rota- 
tion, magnetic activity, accretion or binary interaction that can 
also modulate the observed parameters and evolution of a brown 
dwarf. 



from | B ureasser . McElwain. fe Kirkpatrickl l|2003j) and 
iTinnev et a,U 1120051) : a distance for SDSS 1110+0116 
was estimated using absolute MKO Mj and Mk mag- 
nitude /spectral type relations from iGolimowski et all 
(2004). Dividing the 13 T dwarfs into those with mean 
estimated ages less than 2 Gyr (young) and older than 2 
Gyr (old), we computed the mean ((V ta n)) and standard 
deviation (oy tan ) for each age group. For the young 
sources, (Vtan) = 38 km s _1 and uy tan — 20 km s _1 , 
while for the old sources (Vtan) — 51 km s^ 1 and <rv tan — 
31 km s _1 . The larger mean motion and greater scatter 
in velocities for the latter group is consistent with an 
older mean age. However, a rigorous examination of 
the three dimensional velocities is required before a 
conclusive assessment can be made. 

5.2. Independent Mass and Radius Estimates 

The masses and ages of the objects in our sam- 
ple as derived from the evolutionary models appear to 
be reasonable and consistent with their overall proper- 
ties. However, these values must be considered with 
caution as they are susceptible to systematic errors 
in both the atmosphere and interior models. Dis- 
agreements in mass and radius estimates between these 
two types of models have been suggested in a few 
voung systems IpVlohantv. Javwardhana fc Basril [2004b; 



^lose et al.ll2005l however, see Luhman, Stauffer & Ma- 
majek 2005), and such systematic deviations may be 
present at late ages as well. Fortunately, masses and 
radii can be determined independently of the evolution- 
ary models for those brown dwarfs with bolometric lu- 
minosity measurements, as follows. 

The surface gravity of a solid body, g — GM/R 2 , is 
applicable for brown dwarf photospheres since the verti- 
cal scaleheight of this region (a few km; Griffith & Yelle 
1999) is insignificant compared to the radius of the brown 
dwarf itself (~0.1 R© = 6.95xl0 4 km). Combining this 
with the definition of T e //, L = 4irR 2 o-T^ ff , yields 
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These equations rely only on the T e ff and g values ob- 
tained from the spectral models and the measured lumi- 
nosities, and not on any evolutionary model (c.f., Mo- 
hanty, Jayawardhana & Basri 2004) 

Luminosities for field brown dwarfs have been compiled 
by a number of studies llGeballe et alJl200li iDahn et al 



1 20021 INakaiima et all 120041: ICushing. Ravner. fc Vacca 
2005); here we focus on the r esults of IGolimowsk i et al, 
l|2004|) and lVrba etafl IpOOl . Seven of the T dwarfs in 



our sample have luminosity determinations from these 
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studies; Golimowski et al. also adopt the T e ff and 
log q d eterminations of Gliese 570D from iGeballe et aP 
( 2001). The corresponding masses and radii derived from 
Eqns.|Hland|3 and listed in Table03 generally agree with 
those derived from the evolutionary models. Figure ITU1 
shows a comparison of mass and radius values derived 
from t he Vrba et al. luminosities to the iBurrows et alJ 
l)1997|) theoretical isochrones. With the exception of 
SDSS 1346— 0031 4 , derived values lie between the 1 and 
10 Gyr isochrones, as expected for a Galactic disk sam- 
ple. More importantly, features in the theoretical brown 
dwarf mass-radius relationship are reproduced, including 
the radius minimum of 0.08 R© for larger masses and the 
trend toward larger r adii for lower mass brow n dwarfs 
ijBurrows et al.lll997t IChabrier fc Baraffell2000j) . These 
agreements are promising, and suggest that brown dwarf 
evolutionary tracks are robust at late ages. 

Are these values truly independent of the evolutionary 
models? Not entirely, since the T e // and log// deter- 
minations hinge on the adopted values for Gliese 570D, 
which themselves are partly dependent on the evolu- 
tionary mode l s. T F jf was derived for this source by 
IGeballe et all J200l " using the integrated observed flux 
(over 0.83-2.52 /mi), a bolometric correction determined 
from spectral models (but consistent with more recent 
empirical determinations; Golimowski et al. 2004), the 
measured parallax of the Gliese 570 system, and a ra- 
dius adopted from evolutionary models. In the last case, 
the radii of brown dwarfs at the age of Gliese 570D are 
predicted to be roughly constant, varying by less than 
20% for masses of 0.02-0.07 M©; are determined by 
well-understood interior physics; and have been empir- 
ically tested at the low-mass end with t ransiting extra- 
solar planets (e.g., IBurrows et all (J2004)). Hence, T e ff 
for Gliese 570D can be considered to be empirically ro- 
bust. Its surface gravity, on the other hand, was de- 
rived solely from evolutionary theory. However, given the 
largely constant radii of old brown dwarfs (R « 0.1 R©), 
and assuming a mass in the brown dwarf range (0.02- 
0.07 Mq) yields log// = 4.7 — 5.3 cm s -2 , consistent with 
the adopted model-dependent value. Hence, the adopted 
T e ff and log g for Gliese 570D are only weakly tied to 
evolutionary models, so that the derived parameters for 
other T dwarfs can provide, at mimimum, semiempirical 
tests of these models. 

Two of the objects in our sample, 2MASS 0415-0935 
and 2MASS 0937+2931, are worth addition a l com ment 
as their luminosit ies from iGolimowski et alJ (|2004fl and 
iVrba et aP l)2004t) are significantly discrepant. As a 
result, Eqns. [5] and [7\ yield very different masses and 
radii for these sources. Golimowski et al. determined 
luminosities for individual brown dwarfs by integrating 
the total measured flux over 0.8-5.0 microns (for those 
sources with measured M-band photometry), assuming 
a Rayleigh- Jeans tail for longer wavelengths (with cor- 
rections for molecular absorption between 4 and 15 /mi), 
and using measured parallaxes. Vrba et al. apply a bolo- 
metric correction as a function of spectral type (com- 
puted by Golimowski et al.) to absolute A-band mag- 
nitudes. For 2MASS 0415-0935, the slight differences 

4 The large radius derived for SDSS 1346—0031 could arise if the 
source is an unresolved binary, although systematic effects cannot 
be ruled out. 



between these methods yields a lower luminosity from 
Golimowski et al. (by a factor of 3.5), resulting in a 
similar mass but a much smaller radius (0.083±0.003 
R©) than that derived from the Vrba et al. measure 
(0. 092+0.004 R w ) The former estimate is outside of 
the IBurrows et all l)1997[) model parameter space. Sim- 
ilarly, Golimowski et al. deduce a higher luminosity for 
2MASS 0937+2931 than Vrba et al., and the correspond- 
ing mass (0.118+0.018 M Q ) and radius (0.114+0.008 R ) 
are well outside of the parameter space encompassed 
by the evolutionary models. In both cases, the Vrba 
et al. luminosities yield values consistent with the mod- 
els. This is intriguing, since both 2MASS 0415-0935 and 
2MASS 0937+2931 have measured M-band photometry, 
and the corresponding luminosities from Golimowski et 
al. are expected to be more accurate. These deviations 
may indicate systematic errors in the luminosity deter- 
minations of either Golimowski et al. or Vrba et al., or in 
our T e ff and surface gravity estimates. This is not en- 
tirely unexpected for the apparently metal-poor T dwarf 
2MASS 0937+2931. However, in order to assess the ori- 
gin of these deviations, and whether they actually indi- 
cate problems in the evolutionary models, the number 
and quality of luminosity measurements for low temper- 
ature T dwarfs must clearly be improved. 

6. DISCUSSION 

6.1. The Temperature Scale of Late-type T Dwarfs 

Disentangling the parameters T e ff and g for T dwarf 
spectra enables a more refined examination of the 
T e ///spectral type relation for these objects, a useful 
function for constraining atmospheric properties as 
well as distance estimation. Typically, these relations 
are tied to luminosity measurements and an assumed 
radius, or ra n ge of ra dii for sources with unknown ages 
iDahn et al.l l2002t iTi nney. Burgasser . fc Kirkpatricld 
l2003HGolimowski et al .ll2004HVrha et a,l.ll2004fl . Studies 
have shown that late-type T dwarfs with identical 
spectral types can exhibit significant differences in their 
estimated T e //s. We formally recognize this as the 
additional influence of surface gravity. 

Figure ITU compares the derived T e //s for the sources 
in our sample to their spectral types. Objects with 
low and moderate surface gravities, log// < 5.1, exhibit 
a tight trend of decreasing T e f f wit h increasing spec- 
tral ty pe, largely consistent with the IGolimowski et ail 
l)2004[) T e ff /spectral type relation. Higher surface 
gravity objects, in particular 2MASS 0034+0523 and 
2MASS 0937+2931, have T e// s that are 150-250 K 
cooler for their spectral types. This behavior can be 
understood by the interplay between T e f / and g on the 
major H2O bands, the depths of which determine in part 
T spectral types. T dwarfs with higher surface gravities 
have weaker H 2 bands, and hence earlier spectral types, 
at a given T e f / ; consequently, they would appear to have 
lower T e ffS for a given spectral type. This gravity trend 
is also present when comparing luminosities, as the old- 
est, most massive brown dwarfs (which have the highest 
surface gravities) can have r adii that are 10-15% smaller 
than 1-3 Gyr brown dwarfs ijBurrows et alJll997j) . Cou- 
pled with 10-20% lower T e //s, old brown dwarfs can be 
up to three times fainter than their yo unger field counter- 
parts. This is precisely the deviation IVrba et alJ ( 2004) 
finds in the luminosity of 2MASS 0937+2931 as com- 
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pared to other T6 dwarfs (Golimowski et al. 2004 find a 
somewhat smaller deviation). Hence, luminosity and/or 
T e ff measurements for a consistently classified sample 
could provide a means of segregating young and old sys- 
tems. 

6.2. Improving Temperature and Surface Gravity 
Determinations 

The method outlined here is in some sense a response 
to the current limitations of the spectral models. As the 
models improve in accuracy, direct spectral comparisons 
should eventually be sufficient to determine the physical 
parameters of field brown dwarfs. On the other hand, 
our spectral index comparison method could also be im- 
proved by using additional empirical calibrators such as 
Gliese 570D. Additional calibrators would enable higher 
order corrections to the model ratios, reducing system- 
atic effects. Such empirical constraints can be taken one 
step further: a sufficiently sampled parameter space of 
calibrator sources could enable the determination of T e / / 
and g values independent of the spectral models. To do 
this, several more companion brown dwarfs with inde- 
pendent age and metallicity determinations, and/or bi- 
nary (particularly eclipsing) systems with measured or- 
bital parameters, are required. While Gliese 570D is 
currently the only such calibrator source in the late T 
dwarf regime, three closely-separated late-type T dwarf 
binaries have been identified for which m ass measure- 
ments are feasible ( Bursas ser et al.ll2003d Burgasser et 
al. in prep.); and one T dwarf binary, Epsilon Indi Bab 
(jScholz et al J 120031 fMcCauehrean et al.l|2004|) is also a 
companion to a nearby 0.8-2 Gyr K5 V star. There are 
also several ongoing searches for wide brown dwarf com- 
panions to nearby and young stars (e.g., Chauvin et al. 
2004, 2005; Neuhauser et al. 2005). Identification of sev- 
eral such calibration stars would provide an empirical 
ladder for determining the physical properties of field 
brown dwarfs, and a critical test for both spectral and 
evolutionary models. 

Can we also extend this technique to earlier spectral 
types; e.g., L dwarfs and early- type T dwarfs? In these 
regimes, spectral energy distributions are strongly af- 
fected by photospheric condensates, giving rise to what 
m any consider to be a fo urth "dust" parameter (e . g., /w 
in lAckerman fc Marlevl (I200U) andlMarlev et all ll2002ft: 
T rr inlTsuiil (J2Q05); ao in lBurrows. Sudarskv fc Hubenvl 
(2005)) which may vary with the photospheric gas prop- 
erties or other secondary effects, such as rotation. In 
principle, the method outlined here could be extended 
into the L and early T dwarf regime by incorporating 
this fourth parameter, employing suitable dust cloud 
models and enlarging the sample of empirical calibration 
stars. In practice, this approach may prove more diffi- 
cult as brown dwarf cloud formation remains a poorly 
understood process (e.g., Helling et al. 2004). Neverthe- 
less, T e ff, g and dust content determinations would be 
particularly useful for understanding the unusual tran- 
sition where L and T dwarfs when photospheric dust 
is rap i dly depleted l|Dahn et alJ 120021: iBurgasser et alJ 
l2002aHKnar>r) et al .1120041) . 

6.3. Applications of the Method 

The ability to determine masses and ages for individ- 
ual brown dwarfs is clearly a boon to statistical studies 



of these objects, particularly the Galactic substellar MF 
and formation history. Constraints on these fundamental 
relations have largely been statistical in nature because 
of the d ifficulty in determining masses and ages for field 
sources (iR.eid et alll99HlChabrieil2003lBurgasserl2004t 
lAllen et al.l l2005|^ Using the method described here, 
both distributions can be built up directly from individ- 
ual sources in a sample, assuming that careful consider- 
ation is made of selection effects. While current samples 
of late-type T dwarfs are too small for a robust analysi s 
of this kind fBurgasser et alJ|2002bt iKnapp et al.l f2004'l , 
searches for cold brown dwarfs from wide and deep near 
infrared surveys such a s the UKIRT Infrared Deep Sky 
Survey 5 ijWarre n 2002) may make MF and age distribu- 
tion measurements feasible in the near future. 

In the more immediate term, T e ff and g measure- 
ments, and corresponding mass and age estimates, are 
useful for identifying and characterizing young, low 
mass objects in young star forming regions. Accord- 
ing to current evolutionary models, T dwarfs with ages 
< 10 Myr can have masses of only a few Jupiter masses 
(Mj up ). Several young cluster brown dwarf candidates 
have been identified in rec ent years, including the late- 
type T dwarf S Orionis 70 l | Zapatero Osorio et al.ll2002t 
iMartin fc Zapatero Osoriol feoOS) . Direct comparison of 
spectral data to theoretical models have suggested a very 
low surface gravity for this source, logg ~ 3.5 — 4.0, in- 
dicative of a young, very low mass (~3 WLj up ) brown 
dwarf which has not yet fully contracted. However, such 
direct spectral comparison has been sh own to be flawed 
for lat e T spectral types in general ijBurgasser et alJ 
l2004a|K A calibrated spectral ratio comparison could 
provide a more robust assessment of the physical prop- 
erties of this and other low mass candidates, and verify 
their cluster membership. Furthermore, with luminos- 
ity measurements, independent determinations of mass 
and age would provide semiempirical constraints on the 
evolutionary models from which mass estimates are cur- 
rently derived. Such independent assessments are nec- 
essary to validate the existence of so-called "planetary- 
mass" brown dwarfs in these young star forming regions. 

7. SUMMARY 

We have devised a method for measuring the effective 
temperatures and surface gravities for the lowest lumi- 
nosity T-type brown dwarfs, by the comparison of cal- 
ibrated spectral ratios measured on low resolution near 
infrared spectral data and theoretical models. Using this 
method, we have derived T e ff and g estimates for 14 
T5.5-T8 field brown dwarfs, and a subsolar metallicity 
estimate for the peculiar T dwarf 2MASS 0937+2931. 
Two other sources in our sample, 2MASS 0939-2448 
and 2MASS 1114-2618, appear to have T e// < 700 K, 
and are potentially the coldest brown dwarfs currently 
known. We find no evidence of systematic effects in our 
method, and the agreement of our T e ff determinations 
with prior studies suggests that our results are robust. 
We also find that the scatter observed in T e // /spectral 
type relations likely arises from surface gravity effects, 
as higher surface gravity objects have a lower T e ff at a 
given spectral type. Masses, radii and ages are derived 
for objects in our sample using the IB ur rows et all l)1997fl 

5 See |http:// www.ukidss.org/index.html 
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evolutionary models, while independent mass and radius 
determinations are made for eight T dwarfs with lumi- 
nosity measurements. Broad agreement between these 
values suggests that current brown dwarf evolutionary 
models are accurate at late ages, although this must be 
verified through improved luminosity determinations. 

The comparative technique described here is a useful 
tool for determining the physical properties of the lowest- 
luminosity brown dwarfs, making efficient use of low res- 
olution, and therefore more sensitive, spectroscopy. As 
such, it is a promising method for characterizing large 
samples arising from deep surveys, enabling a direct de- 
termination of the Galactic substcllar MF and formation 
history. While our method remains tied to the current 
generation of spectral models, susceptible to persistent 
opacity deficiencies, the increased use of fiducial calibra- 
tors will ultimately enable a wholly empirical approach, 
allowing critical tests of atmospheric and evolutionary 
theories in addition to the characterization of individual 
brown dwarfs in the vicinity of the Sun. 
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Fig. 1. — SpcX prism spectra for newly observed T dwarfs. All data are normalized at the 1.27 /im flux peaks and offset by a constant 
(dotted lines). Major spectral features are labelled, and regions of strong telluric absorption are indicated by ffi symbols. 
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Fig. 2. — Surface gravity and mctallicity features in T dwarf spectra. The left panel compares the normalized spectra of the T6-T6.5 dwarfs 
2MASS 0937+2931, SDSS 1346-0031 and 2MASS 2228-4310 superimposed on that of the T6 spectral standard SDSS 1624+0029 (red line). 
The right panel compares the normalized spectra of the T7.5-T8 dwarfs 2MASS 0939-2448, 2MASS 1114-2618 and 2MASS 1217-0311 
with that of the T7.5 companion brown dwarf Gliese 570D. Major spectral features are labelled. Note in particular the discrepancies at 
the 1.05 and 2.1 (im peaks due to differences in K I and H2 absorption, respectively. 



14 



Burgasser, Burrows & Kirkpatrick 




Fig. 3. — Comparison of temperature, surface gravity and metallicity effects in T dwarf spectral models. The top panel displays three 
solar metallicity models with logg = 5.0 cm s -2 and T'eff = 1100 (blue dashed line), 900 (black solid line) and 700 K (red dot-dashed 
line). The middle panel displays three solar metallicity models with T e ff = 800 K and logcj = 4.5 (blue dashed line), 5.0 (black solid line) 
and 5.5 cm s — 2 (red dot-dashed line). The bottom panel displays three Teff = 800 K, logcj = 5.5 cm s — 2 models with metallicities [M/H] 
= (black sold line), —0.5 (red dot-dashed line) and —1.0 (blue dashed line). Each spectral model has been deconvolved to the resolution 
of the SpcX data (A/AA ~ 150) and normalized at their J-band peak. 
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Fig. 4. — Comparison of models to spectra in the 1.5—2.4 /im spectral region. Shown are data (solid black lines) for the T6-T6.5 dwarfs 
SDSS 1624+0029, 2MASS 0937+2931 and 2MASS 2228-4310 (left panel) and the T7.5-T8 dwarfs Gliese 570D, 2MASS 0939-2448 and 
2MASS 1217—0311 (right panel). Resampled solar metallicity spectral models for logg = 4.5 (blue dashed line), 5.0 (black solid line) and 
5.5 cm s -2 (red dot-dashed line) and T e ff = 800 K (left panel) and 1000 K (right panel) are overplotted. Both data and models are 
normalized at the H-band peak. 
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Fig. 5. — Values for the spectral ratios H2O— J and K/H as measured on solar metallicity models (points with dash-dot lines connecting 
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Fig. 6.— Derived physical parameter phase spaces for the T dwarfs 2MASS 0034+0523, 2MASS 0050-3322, 2MASS 0243-2453, 
2MASS 0415-0935, 2MASS 0727+1710, 2MASS 0937+2931, 2MASS 0939-2448 and SDSS 1110+0116. T e// and log g values for which 
measurements of the spectral ratios H2O — J and K/H match scaled values for the models (assuming an uncertainty of 10%) are indicated 
by hatched regions. The overlap cross-hatched regions represent our best match "fits" for these parameters. 
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Fig. 7.— Same as Figure|Ufor the T dwarfs 2MASS 1114-2618, 2MASS 1217-0311, 2MASS 1231+0847, SDSS 1346-0031, Gliese 570D, 
SDSS 1624+0029, SDSS 1758+4633 and 2MASS 2228-4310. 



The Physical Properties of T Dwarfs 



19 





Fig. 8. — Teff and logg parameter fit for 2MASS 0937+2931 for spectral models with metallicitics -0.5 < [M/H] < 0, interpolated in 
steps of 0.1 dex. 
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Fig. 9. — Best-fit T e ff and logg values derived for the T dwarfs as compared to evolutionary models from iBurrows 'eTaTI tl997l) . 
Isochrones are indicated by dashed lines (10, 5, 2 and 1 Gyr from top to bottom), while constant mass values (labelled, in Solar masses) are 
denoted by solid lines. Error bars on the data points are based on the breadth of the T e ff, l°g 9 space spanned by the intersecting regions 
in Figures 171 and IS1 (plus an additional 10 K and 0.05 dex uncertainty in T e ff and log g for sampling uncertainty); possible systematic 
errors of 50—100 K in T e ff and 0.1-0.25 dex in logg are not included. The gray circle denotes parameters for 2MASS 0937+2931 assuming 
[M/H] = -0.2 (see § 4.2.2). 
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Fig. 11. — Derived T e ffS for objects in our sample as a function of near infrared spectral type. Sources have been segregated into those 
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and moderate surface gravity objects form a tighter trend than the full sample, as traced by the |Golimowski ct al. 12003) T e / / /spectral 
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ZlvlAoo JU410iy04 — UyoOUOO 


1 55 


U4 10 


1 OS 7/1 

iy . 04 


— uy oo uo. o 


1 7 7n_l_n 
10. ( U±U. 


.uo 


1 7 7/t _i_n 
10.04±U. 


1 1 

11 


17/1 Q_l_n on 
10.4o±U.ZU 


1 7/1 1 ' > 
1 f 4±o 


Z.Z00±U.UUo 


Q A 
0,4 


01\/T A GC 107071 QO/1 1 1 71 OH1 O 

ZMAoo JU 1 z / 15Z4+1 ( 1UU1Z 


T7 


U f z / 


18! 24 


+ 1 / 1U Ul. Z 


lo.6U±U. 


.06 


15.76+0. 


17 


1 7 cfi j_n 1 
lo.o6±U.iy 


110+2 


1 on7J_n nn7 
l.z9/±U.UUo 


3,4 


2MASS J09373487+2931409 


T6p 


09 ,l 37 m 


34! 87 


+29°31'40'.'9 


14.65±0. 


.04 


14.70±0. 


07 


15.27±0.13 


163±4 


1.622±0.007 


3,4 


2MASS J09393548-2448279 


T8 


og^" 1 


35? 48 


- 24° 48' 27'.' g 


15.98±0. 


.11 


15.80±0. 


15 


> 16.6 




1.15±0.06 


2 


SDSS J111010.01+011613.1 


T5.5 


ll fc 10 m 


10!01 


+01°16'13'/0 


16.34+0 


.12 


15.92±0. 


14 


> 15.1 




0.34±0.10 


2,5 


2MASS J11145133-2618235 


T7.5 


11 h 14 m 


51=33 


-26°18'23'/5 


15.86±0 


.08 


15.73±0. 


12 


> 16.1 




3.05+0.04 


2 


2MASS J1217110-0311131 


T7.5 


12 h 17 m 


11-10 


-03° 11' 13'.' 1 


15.86±0 


.06 


15.75±0. 


12 


> 15.9 


91±2 


1.0571±0.0017 


6,7 


2MASS J12314753+0847331 


T5.5 


12 h 31 m 


47! 53 


+08°47'33'.'l 


15.57±0 


.07 


15.31±0. 


11 


15.22±0.20 




1.61±0.07 


1,2,8 


SDSS J134646.45-003150.4 


T6.5 


13 ,l 46 m 


46! 34 


-00°31'50'.'l 


16.00+0 


.10 


15.46±0. 


12 


15.77±0.27 


68±2 


0.516+0.003 


9,7 


Gliese 570D 


T7.5 




14!g6 


-21°21'47'/7 


15.32+0 


.05 


15.27±0. 


09 


15.24+0.16 


169.3±1.7 


2.012+0.004 


10,11 


SDSS J162414.37+002915.6 


T6 


16 h 24 m 


14! 36 


+00°2g'15'.'8 


15.4g±0 


.05 


15.52±0. 


10 


> 15.5 


92+2 


0.3832±0.0019 


12,13 


SDSS J175805.46+463311.9 


T6.5 


17 h 58 m 


05! 45 


+46°33'0g'.'g 


16.15±0 


.09 


16.25±0. 


22 


15.47+0.19 






8 


2MASS J22282889-4310262 


T6 


22 h 28 m 


28! 8g 


-43°10'26'.'2 


15.66±0 


.07 


15.36±0. 


12 


15.30+0.21 




0.31±0.03 


14 



References. — (11 IBurgass er et al.l I 2004bll: (21 ITinnev e t"afl 1200511: (3) IBurgasser et all l2002bH: (4) IVrba eTai] <200l) : (51 IGeballe et all 1200211: (6) 
IBurgasser et all ll93HI): ( 7) iTinn^^lBu rgasser fc Kirkpatrickl 120031) : (8) IKnapaTt^?T72004l) : " iTsvetano^^TTT^ OOOll: f 1 01 IBurgasse^^Ll 1200011 : (11) 
HIPPARCOS IESAI1997I1 : (121 IStrausTet^ril999ll : ( 131 IDahnetail I2002H : (141 IBurgasser. McElv^T7rKirkSrickn200^ . 

a Right Ascension (a) and declination (5) at equinox J2000 from the 2MASS All Sky Data R.clcasc Point Source Catalog ICutri ct al .■2003H . 
^Discovery, parallax and proper motion references. 



TABLE 2 
Log of New SpeX Observations. 



Source 


UT Date 


sec z 


tint (s) 


Calibrator Star 


SpT 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


2MASS 0050-3322 


2004 Sep 07 


1.69 


1440 


HD 225200 


AO V 


2MASS 0727+1710 


2004 Mar 10 


1.00 


1280 


HD 56386 


AO Vn 


2MASS 0937+2931 


2004 Mar 11 


1.02 


720 


HD 89239 


AO V 


2MASS 0939-2448 


2004 Mar 12 


1.43 


1080 


HD 81694 


AO V 


SDSS 1110+0116 


2004 Mar 11 


1.14 


1800 


HD 97585 


AO V 


2MASS 1114-2618 


2004 Mar 12 


1.75 


1080 


hd gsg4g 


AO V 


2MASS 1217-0311 


2004 Mar 11 


1.13 


720 


HD 10g309 


AO V 


SDSS 1346-0031 


2004 Mar 12 


1.09 


720 


HD 116960 


AO V 


SDSS 1624+0029 


2004 Mar 12 


1.06 


720 


HD 136831 


AO V 


SDSS 1758+4633 


2004 Jul 23 


1.22 


720 


HD 158261 


AO V 



Burgasser, Burrows & Kirkpatrick 



TABLE 3 
Spectral Ratios. 



Source 


SpT 


H 2 0- J 


H2O— H 


K/J 


K/H 


Y/J 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


2MASS 0034+0523 


T6.5 


0.103 


0.229 


0.100 


0.218 


0.456 




T7 


0.104 


0.266 


0.180 


0.388 


0.363 




T6 


0.145 


0.297 


0.197 


0.406 


0.444 




T8 


0.041 


0.183 


0.131 


0.255 


0.382 


2MASS 0727+1710 


T7 


0.085 


0.224 


0.164 


0.351 


0.427 


2MASS 0937+2931 


T6p 


0.151 


0.316 


0.076 


0.174 


0.539 


2MASS 0939-2448 


T8 


0.038 


0.149 


0.059 


0.117 


0.493 


SDSS 1110+0116 


T5.5 


0.152 


0.303 


0.217 


0.379 


0.497 


2MASS 1114-2618 


T7.5 


0.039 


0.177 


0.076 


0.150 


0.482 


2MASS 1217-0311 


T7.5 


0.066 


0.207 


0.179 


0.366 


0.374 


2MASS 1231+0847 


T5.5 


0.181 


0.271 


0.157 


0.328 


0.451 


SDSS 1346-0031 


T6.5 


0.131 


0.278 


0.156 


0.351 


0.430 


Gliese 570D 


T7.5 


0.063 


0.198 


0.116 


0.253 


0.397 


SDSS 1624+0029 


T6 


0.154 


0.280 


0.142 


0.311 


0.422 


SDSS 1758+4633 


T6.5 


0.101 


0.247 


0.200 


0.411 


0.400 


2MASS 2228-4310 


T6 


0.157 


0.293 


0.204 


0.440 


0.383 


Corrections' 1 




1.173 


1.567 


0.952 


1.064 


0.883 



a Scale factors applied to model ratios to bring them into agreement with 
measurements for Gliese 570D assuming T^eff = 800 K and logg = 5.1; 
See 8 4.1. 



TABLE 4 

Derived Physical Parameters For T Dwarfs. 



Published Values 



Source 


SpT 




log 9 


[M/H] 






log g c 






(K) 


(cm s~ 2 ) 




(K) 


(K) 


(cm s — 2 ) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


2MASS 0034+0523 


T6.5 


820-860 


5.4-5.5 


-0.2~0 








2MASS 0050-3322 


T7 


960-1000 


4.8-5.0 











2MASS 0243-2453 


T6 


1040-1100 


4.8-5.1 





825-1150 


950-1170 


4.5 


2MASS 0415-0935 


T8 


740-760 


4.9-5.0 





600-750 


690-850 


5.0 


2MASS 0727+1710 


T7 


900-940 


4.8-5.0 





725-950 


830-1020 


4.5 


2MASS 0937+2931 


T6p 


780-840 d 


5.3-5.5 d 


-0.4 0.1 


725-1000 


700-850 


5.5 


2MASS 0939-2448 


T8 


<700 












SDSS 1110+0116 


T5.5 


1020-1100 


4.9-5.2 









4.5 


2MASS 1114-2618 


T7.5 


<700 












2MASS 1217-0311 


T7.5 


860-880 


4.7-4.9 





725-925 


820-1000 


4.5 


2MASS 1231+0847 


T5.5 


1040-1100 


5.2-5.5 









5.0 


SDSS 1346-0031 


T6.5 


960-1020 


5.0-5.2 





875-1200 


950-1180 


4.5 


Gliese 570D 


T7.5 


780-820 


5.1 





784-824° 




5.0-5.3 6 


SDSS 1624+0029 


T6 


980-1040 


5.3-5.4 





800-1100 


920-1100 


5.0 


SDSS 1758+4633 


T6.5 


960-1000 


4.7-4.9 









4.5 


2MASS 2228-4310 


T6 


1080-1140 


4.6-5.0 












a T e f f from Golimowsk i et alJ J20 
b T e// fro mlVrba et alJ PkH). 
c logg from Kn app et alJ 12004 V 
d T e// and logg for [ M/Hl = -0.2. See 8, 4.2.2. 
c T eff and logg from lGeballe et alJ 120011) . 
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TABLE 5 

Estimated Masses, Radii and Ages for T Dwarfs. 



Evolutionary Models 3, Golimowski et al. Luminosities Vrba et al. Luminosities 



Source 


op 1 


M 

jyr 


Lj 
rv 


Age 


l°g l^bol 


M 


L> 

rVj 


log ^bol 


M 


LJ 












IT 

1+0 ) 


\ m 0) 


(.ft©; 


(T 

1+0,1 


\ m @) 


(R^l 

v 41 © ) 




(21 


(31 


(41 


(51 


(61 


(71 


(81 


(91 


(101 


(111 
V- 1 - 1 ) 


2MASS 0034+0523 


T6.5 


0.039-0.055 


0.081-0.090 


3.4-6.9 














2MASS 0050-3322 


T7 


0.022-0.043 


0.090-0.104 


0.5-2.5 














2MASS 0243-2453 


T6 


0.024-0.041 


0.092-0.106 


0.4-1.7 


-5.08+0.06 


0.021+0.002 


0.086+0.005 


-5.03+0.10 


0.023+0.003 


0.091+0 


2MASS 0415-0935 


T8 


0.022-0.044 


0.085-0.101 


1.0-4.9 


-5.73+0.05 


0.020+0.001 


0.083+0.002 


-5.18+0.10 


0.028+0.002 


0.099+0 


2MASS 0727+1710 


T7 


0.022-0.035 


0.093-0.104 


0.5-2.2 


-5.35+0.05 


0.020+0.002 


0.083+0.003 


-5.26+0.10 


0.024+0.002 


0.092+0 


2MASS 0937+2931 


T6p 


0.047-0.063 


0.078-0.084 


5.5-10 


-5.28+0.05 


0.118+0.018 


0.114+0.008 


-5.57+0.08 


0.061+0.009 


0.081+0 


SDSS 1110+0116 


T5.5 


0.028-0.050 


0.087-0.101 


0.5-3.4 














2MASS 1217-0311 


T7.5 


0.022-0.038 


0.091-0.103 


0.7-2.9 


-5.32+0.05 


0.020+0.001 


0.094+0.002 


-5.30+0.09 


0.021+0.001 


0.096+0 


2MASS 1231+0847 


T5.5 


0.038-0.071 


0.078-0.093 


1.6-9.1 














SDSS 1346-0031 


T6.5 


0.031-0.056 


0.082-0.097 


1.0-4.9 


-5.00+0.06 


0.051+0.006 


0.106+0.006 


-5.18+0.12 


0.049+0.006 


0.103+0 


Gliese 570D 


T7.5 


0.041-0.043 


0.087-0.089 


3.7-4.5 


-5.53±0.05 b 


0.037+0.003 


0.090+0.004 








SDSS 1624+0029 


T6 


0.054-0.060 


0.080-0.084 


4.3-5.8 


-5.16+0.05 


0.065+0.007 


0.084+0.005 


-5.11+0.08 


0.073+0.009 


0.089+0 


SDSS 1758+4633 


T6.5 


0.019-0.030 


0.097-0.111 


0.3-0.9 














2MASS 2228-4310 


T6 


0.018-0.034 


0.096-0.115 


0.2-0.9 















a Estimatcs derived from the solar metallicity evolutionary models o f IBurrows et alj 119971) and T e ff and log q determinations from Table 141 
h L bol from lGeballe et all I2001D . 



